Phase transitions in the g-voter model with two types of stochastic driving. 
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In this paper we study nonlinear g-voter model with stochastic driving on a complete graph. We 
investigate two types of stochasticity that, using the language of social sciences, can be interpreted 
as different kinds of nonconformity. From a social point of view, it is very important to distinguish 
between two types nonconformity, so called anti-conformity and independence. A majority of works 
suggests that these social differences may be completely irrelevant in terms of microscopic modeling 
that uses tools of statistical physics and that both types of nonconformity play the role of so called 
'social temperature'. In this paper we clarify the concept of 'social temperature' and show that 
different type of 'noise' may lead to qualitatively different emergent properties. In particularly, we 
show that in the model with anti-conformity the critical value of noise increases with parameter q, 
whereas in the model with independence the critical value of noise decreases with the q. Moreover, 
in the model with anti-conformity the phase transition is continuous for any value of q, whereas in 
the model with independence the transition is continuous for q < 5 and discontinuous for q > 5. 

PACS numbers: 64.60. Ht, 05.70.Ln, 05.50.q 



I. INTRODUCTION 

Recently, various microscopic models of opinion dy- 
namics has been proposed and widely studied by physi- 
cists and social scientists (for reviews see [H-13])- In a 
world of social studies this kind of modeling is known 
as the Agent-Based Modeling (ABM). It has been noted 
recently that despite the power of ABM in modeling 
complex social phenomena, widespread acceptance in the 
highest-level economic and social journals has been slow 
due to the lack of commonly accepted standards of how 
to use ABM rigorously 0, Q . As has been pointed out by 
Macy and Wilier [3| , one of the main problems in the field 
of social simulations is 'little effort to provide analysis of 
how results differ depending on the model designs'. 

The similar problem is visible in a field of sociophysics. 
For example, to study opinion dynamics under confor- 
mity (one of the major paradigm of social response), 
a whole large class of models based on binary opinions 
S = ±1 has beenproposed, among them voter model [q], 
majority ru le M , 3 , Sznajd model Q or nonlinear voter 
models [l(| (HI- For all these models the ferromagnetic 
state is an attractor [l|. On one hand, this is expected 
since the conformity is the only factor influencing opin- 
ion dynamics in these models. On the other hand this is 
obviously not realistic for real social systems. To make 
models of opinion dynamics more realistic several mod- 
ifications has been proposed, among them contrarians 
[H, EH , inflexibles [lij and zealots [l5[ • From the social 
point of view all these modifications describe another ma- 
jor paradigm of social response - so called nonconformity 
1191 . There are two widely recognized types of nonconfor- 
mity: anti-conformity and independence. From a social 
point of view, it is very important to distinguish between 
these two types of nonconformity [l6l Il7j . The term 'in- 
dependence' implying the failure of attempted group in- 
fluence. Independent individuals evaluate situations in- 
dependently of the group norm. From these point of view 



both zealots introduced by Mobilia [l5| , as well as inflexi- 
bles introduced by Galam [3] describe particular type of 
independent behavior. On the contrary, anti-conformists 
are similar to conformers in the sense that both take cog- 
nizance of the group norm - conformers agree with the 
norm, anticonformers disagree. Therefore, contrarians 
introduced by Galam in II 21] or stochastic driving pro- 
posed by Lama et. al. |13| describe anti-conformity. 

Although differences between two types of nonconfor- 
mity are very important for social scientists, the results 
obtained so far indicate that differences may be irrelevant 
from the physical point of view. Both, contrarian and in- 
dependent behaviors, play the role of social temperature 
that induce an order-disorder transition [IH, HH [H, [ll| . 
However, addressing the problem posed by Macy and 
Wilier [3] we would like to check rigorously the differences 
between two types of nonconformity under the frame- 
work of a possibly general model of opinion dynamics. 
In a class of models with binary opinions such a gen- 
eral model has been recently introduced in (TT| under 
the name of g-voter model. As special cases this model 
consists of both linear voter, as well as Sznajd model. In 
this paper we investigate this model in the presence of 
different types of nonconformity and check if results for 
anti-conformity and independence are qualitatively the 
same, according to our first expectation. It should be 
mentioned that another general class of opinion dynam- 
ics, known as a majority rule @,H|, would also be a good 
candidate to test the differences between two types of 
nonconformity. However, introducing a general type of 
independence seem not so straightforward in this case. 

The paper is organized as follows. In the next sec- 
tion we introduce the generalized model with two types 
of nonconformity on a complete graph (topology which 
is particularly convenient for analytical calculations). In 
Section IIIII we analyze the time evolution of the system 
described by the master equation. In this section we 
will already see differences between models with anti- 
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conformity and independence, on contrary to our first 
prediction. In Section IIVI we calculate analytically the 
stationary values of public opinion in a case of infinite 
system. Results presented in Sections IIIII and IIVI indi- 
cate clearly a phase transition between phases with and 
without majority. Therefore in the section [V] we find 
the phase diagram and calculate the transition point as 
a function of model's parameters. Results presented in 
this section show clearly important qualitative differences 
between two types of nonconformity. In the section IVII 
we apply the approach that has been used to study non- 
equilibrium systems with two (Z2) - symmetric absorbing 
states in [I(Jl2l[ to understand more deeply these differ- 
ences and in particularly the origin of the discontinuous 
phase transition in the case of nonconformity . We con- 
clude the paper in the last section. 



II. MODEL 



conformist. This type of anti-conformity has been 
investigated in a case of the Sznajd model on a com- 
plete graph in [19j . Because results depend only 
on the ratio p = pi/pi, in this paper we consider 
Pi = 1 and pi — p. In this version of the model the 
following changes are possible: 
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We consider a set of N individuals on a complete 
graph, which are described by the binary variables 
S = ±1. At each elementary time step q individuals 
Si, . . . , S q (denoted by f for Si = 1 or 4- for Si = — 1, 
where i — 1 , . . . , q) are picked at random and form a 
group of influence, lets call it g-lobby. Then the next 
individual (ff or JJ-), on which the group can influence is 
randomly chosen, we call it voter. 

The part of a model described above is a special case of 
nonlinear q- voter model introduced in (]~Tj | . In the orig- 
inal q- voter model, if all q individuals are in the same 
state, the voter takes their opinion; if they do not have a 
unanimous opinion, still a voter can flip with probability 
e. For q = 2 and e = the model is almost identi- 
cal with Sznajd model on a complete graph [22J. The 



only difference is that in the q- voter repetitions in choos- 
ing neighbors are possible. In [23[ g-voter model with 
e = and without repetition has been considered on a 
one-dimensional lattice. In this paper we also deal with 
a g-voter model with e = and without repetition, but 
additionally we introduce a certain type of noise to the 
model. Original voter model describes only conformity, 
whereas the noise is introduced to describe nonconfor- 
mity. 

In our model conformity and anti-conformity take 
place only if the g-lobby is homogeneous i.e all q individ- 
uals are in the same state. In a case of conformity voter 
takes the same decision as the g-lobby (like in the original 
q- voter model), whereas in a case of anti-conformity the 
opposite opinion to the group. In a case of independent 
behavior, voter does not follow the group, but acts inde- 
pendently - with probability 1 /2 it flips to the opposite 
direction, i.e. S q +± —> —S g +\. 

To check the differences in results that are caused by 
different types of nonconformity we consider three ver- 
sions of the model: 

• Anti-conformity I - with probability p% voter 
behaves like conformist and with pi like anti- 



In other cases nothing changes. 



• Anti-conformity II - with probability p voter be- 
haves like anti-conformist and with 1 — p like con- 
formist. This type of anti-conformity has been in- 
vestigated in a case of the Sznajd model on several 
networks in 11311 and results were qualitatively the 
same as in [19j |. Indeed it is quite easy to notice 
that Anti-conformity II is a special case of Anti- 
conformity I. However, for the record we consider 
here both cases and show that indeed differences 
are only quantitative. In this case the following 
changes are possible: 
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In other cases nothing changes. 



• Independence - with probability p voter behaves 
independently and with 1 — p like conformist. In 
a case of independent behavior individual changes 
to the opposite state with probability 1/2. The 
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following changes are possible: 



Anti-conformity II: 
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• Independence: 



In other cases nothing changes. 



III. TIME EVOLUTION 

In a single time step A t , three events are possible - 
the number of 'up' spins Nf increases or decreases by 1 
or remains constant. Of course all three events can be 
rewritten for the number of 'down' spins N± as Nf+N± = 
A. Also concentration c = N^/N of spins 'up' increases 
or decreases by Ajy = 1/A or remains constant: 



7 +(c) 
7 - (c) 
7 °(c) 



Pro6{c -> c + A N } 
Prob{c — > c - A N } 

Prob {c -)■ c} = 1 — 7 + (c) — 7~(c). (4) 
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Whereas for the infinite system the above formulas 
take much simpler forms: 

• Anti-conformity I: 

7+ = (1 - c)c q + p(l - c) 9+1 

7 - = C (l - c) 3 + P C 9+1 



(9) 



The time evolution of the probability density function 
of c is given by the master equation: 

p(c,t + A t ) = 7 + (c — Ajy)p(c — Ajy, t) 
+ j-(c + A N )p(c + A N ,t) 
+ [1-7 + (c)-7-(c)]pM). (5) 

Of course analogous formula can be written for Nf. The 
exact forms of probabilities 7 + (c) = 7 + (A|) = 7+ and 
7~(c) = 7~(Af) = 7 _ depend on the version of the 
model and for the finite system are the following: 

• Anti-conformity I: 
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• Anti-conformity II: 

7+ = (l-p){l-c)c q +p(l-c) q+1 
7" = (1 - p)c(l - c) q + pc q+1 

• Independence: 

7 + = (l-p)(l-c)c 9 +p(l-c)/2 
7" = (1 -p)c(l - c) q +pc/2 



(10) 



(11) 



Solving analytically master equation ([5]) is not an easy 
task, but exact formulas for 7+ and 7" allow for a nu- 
merical solution of the equation. For an arbitrary initial 
state the system reaches the same steady state. In a case 
of the infinite system the probability density function is 
a sum of delta functions p s t(c) — S(c — Ci) + S(c — C2) + 
. . . + 8(c — Cfc), whereas for the finite system p s t(c) has 
maxima (peaks) for the c = Cj j = 1, . . . , k which are 
getting higher and more narrow with the system size, 
approaching deltas for the infinite system. The number 
of peaks k and values c\ , . . . , Ck depend on the version of 
the model, as well as on the model's parameters p and q. 

Examples of the stationary probability density func- 
tions for the g-lobby of size q — 7 and the system size 
A = 200 are presented in Figs. [T] (anti-conformity I) 
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FIG. 1: Stationary probability density function of the concen- 
tration of 'up' spins for the q- voter model with anti-conformity 
I for the system of N = 200 individuals and the size of the 
lobby q — 7. As seen for small values of anti-conformity 
p the system is polarized, but for large values of p there is 
no majority in the system. For p = (the case without 
anti-conformity) the system consists of all spins 'up' or all 
spins 'down'. With increasing p maxima are getting lower 
and approaching each other. Eventually they form a single 
maximum. This is typical behavior for a continuous phase 
transition. The critical value of p can be found analytically 
(see section |V)) and depends on q. 



FIG. 2: Stationary probability density function of the concen- 
tration of 'up' spins for the g-voter model with independence 
for the system of TV = 200 individuals and the size of the lobby 
q = 7. As seen for small values of independence p the system 
is polarized, but for large values of p there is no majority in 
the system. For p — (the case without independence) the 
system consists of all spins 'up' or all spins 'down'. For larger 
values of p the third maximum appears at c = 1/2 (no ma- 
jority) . This maximum increases with p while the remaining 
two maxima are decreasing. Above certain value of p there is 
only one maximum for c = 1/2. This is typical behavior for a 
discontinuous phase transition for which we can observe the 
phase coexistence. 



and [2] (nonconformity) . As seen from Figs. [T] and [3] for 
small values of noise p (whether the noise is introduced 
as independence or anti-conformity) the system is polar- 
ized, whereas for large values of p there is no majority in 
the system. However, the transition from the phase with 
majority to the phase without majority is very different 
for each type of noise. In a case with anti-conformity we 
observe continuous phase transition for arbitrary value of 
q, whereas in a case with nonconformity there is a contin- 
uous phase transition for q < 5 and discontinuous phase 
transition for q > 5. 

In the case with anti-conformity two states with ma- 
jority, represented by two equally high peaks, are stable 



below the critical value of p* . As p < p* increases the 
two peaks approach each other and eventually for p = p* 
they form a single peak, which is a typical picture for a 
continuous phase transition (see Fig. [IJ HE] ■ In a 
case with independence this picture is valid only for the 
lobby q < 5. For q > 5 the transition is very different. 
Again for small values of p there are two peaks but with 
increasing p they are not approaching each other. In- 
stead, for p = p* the third peak appears at c = 1/2 (see 
Fig. [2]). The third peak is initially lower than remaining 
two peaks, which means that it represents a metastable 
state. As p > pi increases the third peak grows and for 
p = p\ all three peaks have the same hight. For p > p\ 
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FIG. 3: Time evolution of the concentration of 'up' spins for 
the model with anti-conformity with lobby q = 7 and level of 
anti-conformity p — 0.5. The system size N = 200. Sponta- 
neous transitions between two stationary states are visible. 



the central peak dominates over the other two, which 
means that the state c = 1/2 is stable and remaining 
two are metastable. Finally, for p = p^ side peaks disap- 
pear and only the center peak remains. This is a typical 
picture for a discontinuous phase transition, which takes 
place at p = p% [H, [25[. Two values of independence 
parameter p = p* and p = p%, demarcate the existence of 
metastability (spinodal lines) [26|, |28[ . Values p*, p|j Pt 
depend on the size of the lobby q, which will be shown 
exactly in the section Ivl 

Before moving on to the analytical results for the in- 
finite system and determine the points of phase transi- 
tions, let us present the time evolution. We stop for a 
moment on the case of a finite system. Having exact 
formulas for transition probabilities 7 + and 7^ we are 
able not only to calculate numerically stationary density 
function p s t{c) but also generate sample trajectories of 
concentration (Figs|3J |H [3] and H]). In a case of a finite 
system spontaneous transitions between states are possi- 
ble. In a case with anti-conformity transitions between 
two states, which correspond to peaks in the probability 
density function p s t(c), are possible below critical value 
of p. Because both peaks are equally high the system 
spends the same time on average in each state. 

This is also true in the case with independence and q < 
5 (see Fig. [4]). However, as we have already written for 
q > 5 there is a discontinuous phase transition between 
state with and without majority and forp € (pl,p^) there 
are three possible states. Therefore for q > 5 we expect 
spontaneous transitions between three states. 

Such transitions are indeed observed. For p G (p*,pl) 
state with majority is stable and state without majority 
is metastable. Therefore, the system spends more time 
in states with majority. For p € iPziPs) the situation 
is exactly the opposite - the state without majority is 
stable. In a transition point p — p^ all three states are 
stable and the system spends the same time on average 
in each of three state (see Figs. [5] and [5]). 




time(MCS) x1Q 4 

FIG. 4: Time evolution of the concentration of 'up' spins for 
the model with independence with lobby q = 5 and level of 
anti-conformity p — 0.175. The system size TV = 200. Spon- 
taneous transitions between two stationary states are visible. 




FIG. 5: Time evolution of the concentration of 'up' spins for 
the model with independence with lobby q = 7 and level of 
anti-conformity p = 0.111 (for this value all three states are 
stable). The system size TV = 200. Spontaneous transitions 
between three stationary states are visible. 



IV. STATIONARY CONCENTRATION 

In the stationary state we expect that the probability 
of growth 7+ should be equal to the probability of loss 
7~ and therefore: 

F{c,q,p) =7 + (c,(?,p) -j~(c,q,p) = 0, (12) 

where F(c, q 7 p) can be treated as an effective force - 7+ 
drives the system to the state 'spins up', while 7~ to 
'spins down'. Therefore we can easily calculate also an 
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time(MCS) x1Q 5 

FIG. 6: Time evolution of the concentration of 'up' spins 
for the model with independence with lobby q — 9 and level 
of anti-conformity p = 0.0685. The system size TV = 200. 
Spontaneous transitions between three stationary states are 
visible. 
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effective potential: 

V(c,q,p) = - / F(c,q,p)dc. 



Independence 



(13) 



To calculate stationary values of concentration we simply 
solve the equation: 



F(c,g,p)=0, 



(14) 



or alternatively find minima of the potential V. Although 
the first possibility is more straightforward, we will see in 
the next section that knowing the form of potential will 
help us to calculate the transition points. 

Exact forms of the force F and the potential V for the 
infinite system are given below: 

• Anticonformity I 

F = (1 - c)c q + p(l - c) q+1 - c(l - c) q - pc q+1 
V = — (c q+1 + (1 - c) q+1 ) 



+ ^ (c 9+2 + (1 - c) 9+2 ) 

• Anticonformity II 

F = (1 - p)(l - c)c q + p(l - c) q+1 

- (1 -p)c(l -C)« - P C q+1 



(15) 



V 



(c q+1 + (1 - c) q+1 ) 



(16) 



• Independence 

F = (1 - p)(l - c)c q + 



V 



(l-p)c(l-c)«-| 

-1Z£ + (1- C )«+1) 

+ ^pf ( C 9+2 + (1 - C ) 9+2 ) - - c) (17) 



Solving analytically Eq. (IT41) . i.e. finding c s t as a func- 
tion of p for arbitrary value of q is impossible, but we can 
easily derive the opposite relations satisfying Eq. (fT4|) : 

• Anticonformity I 



P = 



c st (l - c st ) q - (1 - c st )c q s 



(18) 



Anticonformity II 



c«t(l ~ c si ) 9 - (1 - c st )c!? 



(1 - c st ) 9+1 + c st (l - c s4 )« - (1 - c st )c 



p = 



C st (l - Cst) 9 - (1 - C st )c q s 



(1 - c st )/2 + c,t(l - c st ) 9 ~ (1 ~ c 5t )c« t - c st /2 



<20) 



We have used the above formulas to plot the depen- 
dence between steady value of concentration c s t and the 
level of noise p for several values of q (see Fig. [7]). Al- 
though, only relation p(c st ) is calculated analytically and 
the opposite relation is unknown, we plot c s t{p) simply 
rotating the figure with the relation p(c st ). The clear dif- 
ferences between two types of noise are visible - in a case 
with anti-conformity the transition value of p increases 
with q and in a case with independence it decreases 
with p. Moreover, the type of transition is the same 
for arbitrary value of q in a case with anti-conformity, 
whereas in a case with independence the transition be- 
tween phases with and without majority changes its char- 
acter for q > 5. 

It should be also noticed that formulas (|18j) - (j20)) has 
been obtained from the condition (1141) . i.e correspond to 
extreme values of potentials (TT5l) - (fTT|) . However, only 
minima of the potential correspond to the stable value 
of concentration. Therefore, in Figs. [7] and [8] we have 
denoted unstable values that correspond to the maxima 
of potentials by the dotted lines. Moreover, we have pre- 
sented the flow diagram for chosen values of q to show 
precisely which state is reached from a given initial con- 
ditions. Particularly interesting behavior is related with 
independence (bottom panel in Fig. [5]). Starting from 
two different initial concentrations disorder or order can 
be reached as a steady state (hysteresis). 

In the next section we derive analytically transition 
points using the knowledge of the effective potentials V 
(US]), (HU) and (HZD. 



V. PHASE TRANSITIONS 

As already noticed there is a continuous phase transi- 
tion for the model with anti-conformity I and II for arbi- 
trary value of q. Below critical value p = p*(q) an effec- 
tive potential has two minima and above only one. Con- 
sequently stationary probability density function p s t{c) 
for p < p* has two maxima and for p > p* only one at 
c = 1/2 (there is no majority in the system). Analogous 
behavior is observed for the model with nonconformity 
but only for q < 5. In all these cases we can easily cal- 
culate critical value p* making the simple observation 
concerning the behavior of the effective potentials (|15p . 
([T6]l and (HU) for q < 5 at c = 1/2 (see also Fig. Q] for 
clarity): 

• For p < p* potentials V(c,p, q) have the maximum 
values for c = 1/2 and therefore 



d 2 V(c, P ,q) 



dc 2 



< 



(21) 
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FIG. 7: Phase diagram for the model with anti-conformity 
(top panel) and independence (bottom panel). Dependencies 
between steady values of concentration c 3 t and the level of 
noise p for several values of q are plotted using formulas (|18|l - 
(I20|l . Although, only relation p(c s t) is calculated analytically 
and the opposite relation is unknown, we plot c a t(p) simply 
rotating the figure. Dotted lines have been used to mark 
instability. Although both types of line (solid and dotted) are 
obtained from equation <|14[) . i.e correspond to extreme values 
of potentials (|15I17[) . only solid lines denotes stable values, i.e. 
correspond to the minima of potentials (see also Fig. [S]). The 
clear difference between two types of noise is visible - in a 
case with anti-conformity the transition value of p increases 
with q and in a case with independence it decreases with p. 
Moreover, the type of transition is the same for arbitrary value 
of q in a case with anti-conformity, whereas in a case with 
independence the transition between phases with and without 
majority changes its character for q > 5. 



FIG. 8: Flow diagrams for the model with anti-conformity for 
q = 2 (top panel) and independence for q = 9 (bottom panel). 
Particular values of q have been chosen just as examples and 
the dependences between stationary values of c and parame- 
ter p for other values of q are seen in Fig. [7] Here solid lines 
denote stable (attracting) steady values of concentration that 
correspond to the minima of potentials (|15l) - (|17[) . whereas 
dotted lines denote unstable values of c that correspond to 
maxima of potentials. Arrows denote the direction of flow, 
i.e. how the concentration changes in time. Particularly inter- 
esting behavior is related with independence (bottom) panel. 
Starting from two different initial concentrations disorder or 
order can be reached as a steady state (hysteresis). 



This means that for p = p* maximum changes to mini- 
mum at c = 1/2: 



For p > p* potentials V(c,p, q) have the minimum 
values for c = 1/2 and therefore 



d 2 V(c,p,q) 







dc 2 

Hence, the critical values: 
• For anti-conformity I: 



dF(c,p, q) 



dc 



0. (23) 



d 2 V{c,p,q) 



> 0. 



(22) 



(24) 
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FIG. 9: Schematic plot of a potential for the model with 
independence and q > 5. For p G (0,p*), potential V(c,p, q) 
has 2 minima that correspond to the states with majority. 
For p G (pijPZ) the potential V(c,p, q) has 3 minima and 
the state without majority is metastable. For p G (p2,P3) the 
potential V(c,p, q) has 3 minima and the states with majority 
are metastable. Finally, for p G (p^, 1) the potential V(c,p, q) 
has only 1 minimum that corresponds to the state without 
majority. Exact form of the potential is given by Eq. (|17[1 . 



• For anti-conformity II: 

p*(q) = g ^ L ( 2 5) 

2q 

• For independence with q < 5: 

p*(q) = q ~\ : (26) 

As we see simple calculations allowed us to find critical 
points for almost all cases, except of the model with non- 
conformity for q > 6. In all cases considered above, there 
is a continuous phase transition between phases with and 
without majority. However, for model with independence 
and q > 6 the phase transition becomes discontinuous 
which has been already discussed in the section Hill This 
behavior can be also suspected from the form of an ef- 
fective potential (|17[) which for q > 6 has the following 
properties (see also Fig. ^j: 

• for p G (0,Pi), V(c,p, q) has 2 minima 

• for p — pi in V(c,p, q) 3-rd minimum emerges 

• for p G (ptiPz)) V(. c ,Pi<l) has 3 minima 

• for p = pi,, V(c,p, q) has 3 minima which are equal 

• for p G (PziPs), V( c ,P,q) has 3 minima 



• for p — in V(c,p, q) 3-rd minimum disappears 

• for p G (pjj, 1), V(c,p, q) has 1 minimum 

As we see, there is an interval p G (p*,p-}) in which po- 
tential V(c,p, q) has 3 minima and therefore stationary 
probability density function has 3 maxima (see Fig. [2]). 
In this region we have a coexistence of two phases - with 
and without majority. For p < pi, the state with major- 
ity is stable and the state without majority is metastable 
and for p > pi, the state with majority is metastable and 
without majority stable. Consequently the phase tran- 
sition appears at p = p\ — p* and p — p\,p% designate 
spinodal lines [26l |28|. 

Transition points p* as a functions of q for all three 
models are presented in Figure 1101 As seen, both mod- 
els with anti-conformity (I and II) behave qualitatively 
the same - critical value of p increases with q. However, 
for the model with independence the transition point de- 
creases with p. Moreover, for q = 5 in a case with inde- 
pendence, phase transition changes its type from continu- 
ous to discontinuous. To clarify our results we decided to 
present the complete phase diagrams for the models with 
anti-conformity and independence in Fig[TT] Because re- 
sults for both models with anti-conformity (I and II) are 
qualitatively the same we present the phase diagram only 
for the model with anti-conformity II. 

The first difference between models with anti- 
conformity and independence connected with qualitative 
dependence between p* and q is easy to explain heuristi- 
cally. This is quite obvious why in the model with inde- 
pendence the critical point p* decreases when q increases. 
When q increases it becomes unlikely to choose randomly 
q parallel spins and therefore the noise term dominates 
because it is independent on the state of the g-lobby. 
Similarly, it can be understood why in the model with 
anti-conformity the critical point p* increases with q. 
It should be recall here that anti-conformity takes place 
only when q+l parallel spins are chosen randomly, which 
is more unlikely than choosing q parallel spins. There- 
fore the anti-conformity term decline in importance even 
more than conformity term as q increases. The second 
difference between models that corresponds to the change 
of the transition type in the model with independence is 
not so easy to understand intuitively. This result has 
been obtain numerically from the potential (|17|) . but in 
the next section we will show how this result could be 
also derived from an approximate Landau description. 

VI. LANDAU DESCRIPTION 

Although we were able to calculate critical points for 
the model with anti-conformity and for the model with 
independence with q < 5 directly from the potentials 
(fT5j) , (IT6l) and (fT7f , it can be instructive to use the classi- 
cal description proposed by Landau for equilibrium phase 
transitions [25]. It has been shown, that this kind of de- 
scription can be also obtained as a mean field approach 
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— □ — anti-conformity I 
0-9 ~ — A — anti-conformity II 




FIG. 10: Transition points p* as a functions of q for all three 
models. Solid lines denote the line of the phase transition 
and dashed lines denotes spinodal lines i.e. determine the re- 
gion with metastability. Several differences between models 
are visible. As seen both models with anti-conformity be- 
have qualitatively the same - critical value of p increases with 
q. However, for the model with independence the transition 
point decreases with p. Moreover, for q > 6 phase transition 
changes its type from continuous to discontinuous. 



for the Langevin equation of non-equilibrium systems 
with two (Z2) - symmetric absorbing states [20I. l2lj. 

In our paper we have written the master equation as 
a function of concentration c = Nf/N of 'up'-spins. We 
have decided to use this quantity for convenience - cal- 
culations are simple and equations have compact forms. 
However, to meet the request of the symmetry [2(| [2l[ 
one should use an order parameter (in this case magne- 
tization) defined as: 



N 



(27) 



for which potentials (fl5|) . (fT6| and (fTTj) are symmetric 
under reversal <f> — ¥ — <j). 

Following the approach presented in [2^, [2l| , which co- 
incides with the classical approach proposed by Landau, 
we expand potentials (fl~5j) . (|T6|) and §17\ . rewritten as 
a function of </), into power series and keep only several 
terms of an expansion: 



V{4>) = A<? 



C4>\ 



(28) 



where coefficients A = A(p,q),B — B(p,q),C = C(p,q) 
depend on the model. 

For the model with independence: 



A(p,q) = - 
B(p,q) = - 
C{p,q) = - 



(1 ~P) (g - 1) P 

2i 2 4' 
(l-p)g(g-l)(g-5) 

2« 24 
(l-p)q(q-l)(q-2)(q-3)(q-9) 



2" 



720 



(.29) 




"a. 0.2 




FIG. 11: Phase diagrams for the model with anti-confomity 
II (top panel) and independence (bottom panel). As seen for 
the model with anti-conformity critical value p* increases and 
for the model with independence decreases with q. For anti- 
conformity (top panel) there is a continuous phase transition 
(denoted by the solid line) between order (i.e. c 7^ 1/2 or 
equivalently m 7^ 0) and disorder (i.e. c = 1/2 or equivalently 
m = 0). In the model with independence (bottom panel) 
there is a continuous phase transition only for q < 5. At 
q — 5 the phase transition changes its type from continuous 
to discontinuous. For q > 5 an area in which one of two 
phases (ordered or disordered) is metastable is limited by so 
called spinodal lines. This area is signed as 'coexistence', 
although the real coexistence occurs only on the transition 
line. However, in the region of metastability both phases can 
be observed depending on the initial conditions (hysteresis), 
which can be also seen from the flow diagram in Fig(8] 



From the Landau theory it is known that for B{p 1 q) > 
and C(p, q) > there is a critical point at which A(p, q) 
changes sign (25[. For A < potential V{4>) has two 
symmetric minima and thus the system is driven to one 
of partially ordered state with 7^ 0. For A > potential 
V((j)) has minimum at <j) — and therefore the system 
remains in an active disordered state and a magnetization 
4> fluctuates around zero. From (|29|) it is easy to calculate 
that: 



A(p, q)=0 -> p = p* 



ff-1 



A < 
A > 



P < P 
P>P* 



) = 0, 



2 g-i 



(30) 



which coincides with the result (|26|) obtained from the 
not an approximate version of the potential (|17[) . 
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As shown within the classical Landau theory, for 
B(p,q) < and C(p,q) > a discontinuous jump in 
the order parameter is expected (25[. Again from ([29| it 
is easy to see that B(p,q) < for q > 5 (see also figure 
I12p . Therefore we expect discontinuous phase transition 
for q > 5, which also agrees with results obtained from 
(fTT|) . It should be mentioned here that transition for 
B < could possibly be included in the class of gener- 
alized voter model (so called unique GV transition) [2(| ■ 
It has been noticed for a general class of models with 
two (Z 2 ) - symmetric absorbing states that for B < 
the location of the potential minimum changes abruptly 
from cf> = to <f> ± 1, i.e. discontinuous phase transition 
is observed [2(|. In our model the situation is slightly 
different, because for p > there are no absorbing states 
and below the transition point \4>\ < 1. However, still 
the system jumps from a totally disordered to a partially 
ordered state, i.e. discontinuous phase transition is ob- 
served. 

One should also notice that in the case of a q-voter 
model with independence for q > 9 also C(p, q) becomes 
negative and than the approximation (|28[) is not valid 
anymore. 

Analogous calculations can be done for the models 
with anti-conformity. Because both models with anti- 
conformity are qualitatively the same we present here re- 
sults for the model with anti-conformity II. In this case: 



A(p, q) 
B(p,q) 



2pq-q + l 
1 1-p 

ip_(q + i 

42i \ 3 



C(p,q) = - 



11 -p 



q-1 



1 



q-i 

3 



6 2« 

lp_(q + 1 

6 29 ^ 5 

Therefore in the case with anti-conformity: 

q-i 



(31) 



A(p,q)=0 

A < 
A > 



P = P = ,, 
2q 

p <p*, <P o 
p>p\ = 0, 



(32) 



which coincides with the result (|25|. Moreover, for the 
model with anti-conformity (see figure I12[) coefficient 
B(p = p* ,q) > for any value of q, i.e. the transition is 
continuous for arbitrary value of q. 



VII. CONCLUSIONS 

In this paper we have asked the questions about the im- 
portance of the type of nonconformity (anti-conformity 
and independence) that is often introduced in models of 
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FIG. 12: Coefficient B(p,q) (see an effective potential (J28J) ) 
for a critical point p = p* at which A(p, q) changes sign. 
For the model with independence (denoted by V) coefficient 
B < for q > 5 which suggests discontinuous phase transi- 
tion, whereas for the model with anticonformity (in the figure 
the model with anticonformity II is presented by '*') coeffi- 
cient B > for any value of q and therefore the transition is 
continuous for arbitrary value of q. 



opinion dynamics (see e.g. [H GJ, M, Q ) ■ We realized 
that the differences between different types of nonconfor- 
mity are very important from social point of view [16j but 
we have expected that they may be irrelevant in terms of 
microscopic models of opinion dynamics. To check our 
expectations we have decided to investigate nonlinear q- 
voter model on a complete graph, which has been recently 
introduced as a general model of opinion dynamics 
To our surprise results for the model with 
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our surprise results tor the model witii anti- 
conformity are qualitatively different that for the model 
with nonconformity. In the first case there is a continu- 
ous order-disorder phase transition induced by the level 
of anti-conformity p. Critical value of p grows with the 
size of the q-lobby. On the other hand for the model with 
nonconformity the value of the transition point p* decays 
with q. Moreover, the phase transition in this case is con- 
tinuous only for q < 5. For larger values of q there is a 
discontinuous phase transition - and coexistence of or- 
dered (with majority) and disordered (without majority) 
phase is possible. 

We have suggested in the title and the introduction 
of the paper that both types of nonconformity play role 
of the noise. However, only independence introduces the 
real random noise that plays the similar role as a temper- 
ature. From this point of view the change of the type of 
transition reminds of the similar phenomena in the Potts 
model (for review see [29j]). In the Potts model there is 
a first-order phase transition for q > 4 and second-order 
phase transition for smaller values of q, where q denotes 
the number of states of the spin. Of course in a case of 
our model q does not denote the number of states, which 
is always 2, but the size of the group. In a way similar 
observation has been recently done by Araujo et. al. [3(| 
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within the model of tactical voting. They have consid- 
ered q candidates on which citizens vote and proposed a 
balance function to quantify the degree of indecision in 
the society due to the coexistence of different opinions. 
It turned out that for some values of model parameters 
the model boiled down to the g-state Potts model, al- 
though similarly like in our model q denoted the number 
of candidates instead of the number of states. The similar 
change of the type of transition has been also observed 
in a general class of systems with two (Z 2 ) - symmetric 
absorbing states within a Langevin description [2(J HH ■ 
Moreover, it has been suggested that models with many 
intermediate states (i.e. the Potts model or a simple 3- 
state model described in [HI) behave as equivalent two- 
state models with effective transitions that are nonlinear 
in the local densities [21], which is the case of q- voter 
model or a two-state model of competition between two 
languages [3l| . The theory presented in [2(], HH suggests 
also that the continuous phase transition that is observed 



for q < 5 could possibly be included in the Ising class, 
whereas discontinuous that is observed for q > 5 in the 
class of generalized voter model. 

Concluding the paper we would like to pay attention 
to one more phenomena that is visible in Fig. [TO] For 
lobby q — 2 results are the same for anti-conformity and 
independence. Therefore there is no surprise that the 
difference between two types of nonconformity has not 
been noticed while studding the Sznajd model (i.e.q = 2) 

mm. 
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